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Abstract 

The decays ij/rj' -^ tttttt and rj' -^ r]Tnr are studied up to leading and next-to-leading order 
within the framework of C^l(3) ® Ur{3) Chiral Perturbation Theory. The analysis incorporates 
important features of the tj — tj' system, such as the contribution of the glueball aGG due to 
the axial anomaly and %/% mixing. One-loop corrections, which are third-order contributions 
according to the combined chiral and 'i/Nc expansion, are not included. Reasonably good results 
are obtained in most cases. 

PACS: 12.39 Fe, 13.25 -k, 14.40 Aq, 14.40 Cs. 
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1 Introduction 

Chiral perturbation theory has proved to be a good tool to describe the dynamics of the low-energy 
region of QCD, where the natural degrees of freedom are the eight Goldstone bosons (pions, kaons 
and 7]) associated with the spontaneous breaking of SUfi{3) ® SUl{3) chiral symmetry. If a large 
number of colors A^'^ is allowed, the theory can be enlarged to the so-called Ufi{3) (X" Ul{3) Chiral 
Perturbation Theory, that includes a ninth particle — the r]' . 

The r] — > vrvrvr decays have drawn a lot of attention from theoreticians since they provide a 
measure of the isospin violation due to the quark masses: the main contribution to the amplitude 
is proportional to md — rriu- In principle, this process can be analyzed in the simpler octet theory. 
However, the rj particle is actually a superposition of the SU{3) singlet r/g and the U{1) singlet %, 
and mixing is seemingly important [9 ~ —20°), so a model including this effect is expected to give 
better results. 

The measured values Q for i] decays are the following: 

r(r? -^ vroTTovro) = 379 ± 40 eV , 
r(r/ ^ 7ro7r+7r_) = 274 ± 33 eV , 



T^^P{rj — > 7ro7r-|-7r_^ 



1.35 ±0.05. 



The U{3) theory can be also used to study rj' decays. Two different channels will be analyzed 
in this paper: r]' -^ vrvrvr and rj' -^ r/vrvr. 

The T]' — > vrvrvr transition is also an isospin violating process. It is however not one of the 
dominant decay channels for the rj' , as happened in the corresponding ry decay. As a consequence, 
the branching ratios associated to these decays are more difficult to determine experimentally 
because they stem from a small fraction of the total observed events, so the uncertainties are 
unfortunately higher: 



r^^p(V - 


-^ vTovroVTo) = 311 ± 77 eV . 


r^^p(r?' - 


> vrovr+vr_) < 1005 eV , 




r > 0.3 



The r]' -^ tjittt transition is in contrast one of the most important decay channels for the r]' . 
The measured rates are: 

r^^P[r]' -^ 7/^ovro) = 42.0 ± 6.0 keV , 
r^^P[r]' -^ rj7T+TT-) = 88.9 ± 10.0 keV , 
r(V -^ ??vrovro) 



T{r]' -^ r/vr+vr_) 



2.1 ±0.5 . 



In general, the description of these decays involves the estimation of the following amplitudes: 

{ViV2V3\Vo)=i(.2Tr)W'{pf-pi)Av,^'P^r2r,{s,t,u) , 



where 

s = (po — Pi) ; t = (po — P2) ', u = {po — P3) {Pi is the momentum of particle Vi) . 

The decay rates in the center-of-mass reference frame require a phase space integral of the squared 
amplitudes over some region TZ, defined by the kinematic restrictions of a three body decay. 

Some general issues can be inferred from symmetry considerations. In the isospin limit (and 
in the absence of electromagnetic interactions) the r]/r]' — > tttttt decays are strictly forbidden by 
Bose symmetry. Therefore, they must originate in isospin-violating terms, which are proportional 
to ruu — m^. One can see Q that the pions must emerge in an / = 1 configuration. This A/ = 1 
selection rule can be used to prove a relation between the two decay channels: 

Aooo(s,t,M) = AQ+_{s,t,u) + AQ+^{t,u,s) + AQ+_{u,s,t) . (1) 

If one assumes that the dominant decay channel is the vectorial one (by vector-meson dominance 
y]) where rj/r]' -^ irp ^ tttttt, one can assume the amplitudes to be flat, moment independent 

functions: Aqh {s,t,u) ss Aq^ (s) « Aq^ (-^p) = A. Under this assumption, the ratio between 

the charged and the neutral channels can be easily estimated: 



T{ri -^ TToTToTTo) _ ^ /-/^ |3^P 



r(r/^7ro7r+7r_) /t^ |^P 



1.5 . (2) 



Notice that phase space corrections due to M^o / M^+ have been neglected (the region of integra- 
tion TZ is the same in the numerator and the denominator). They can easily be included within 
the present approximation, but turn out to increase rather than decrease the ratio — the main 
corrections to it must come from the inclusion of other decay channels. 

The isospin symmetry constraints on the r]' — > tttttt system produce more restrictive results. In 
this case, the wave function must be symmetric under the exchange of pions only. Besides, the 
total isospin for the three-particle state must be equal to zero. When this is taken into account, 
the expansion in terms of Clebsch-Gordan coefficients leads to a simple relation: 

The amplitudes being equal, the rates will be identical except for the combinatorial 1/2! pre- 
factor in the neutral case, so the ratio must be equal to 2. Phase-space corrections produce a 
deviation from this value. 

A first approach to the r/ decays can be done in the usual SU{'i) Chiral Perturbation Theory. 
The leading-order results are definitely too small. The 0{j>^) result is closer to the experimental 
data, because the unitary corrections due to the final-state interactions are surprisingly large, 
especially in the 1=0 and 1, S-wave vrvr channels. These results (collected in section 5) seem to 
indicate that the rj -^ tttttt decays are dominated by the (well-known) vectorial resonance and by 
an intermediate low-mass scalar resonance, the celebrated a particle ^ |[. 

Nevertheless, as mentioned before, there is another possible reason why the SU{3) Chiral Per- 
turbation Theory prediction fails: the physical 77 particle is not one of the states in the octet, but 
a superposition of r/g and rjQ. The rjo, on the other hand, is a mixture of the pseudoscalar quark 
current and the gluonic state GG due to the axial anomaly. All this should be taken into account 
in order to get a more careful description of the ry meson Q] . 



2 Ul{S) (g) Ur{3) Chiral Perturbation Theory 

f^L(3) (S) t/ij(3) Chiral Perturbation Theory is an appropriate tool to deal with the r] — rj' — it 
transitions. In the chiral limit ruq = 0, the QCD Lagrangian is invariant under the flavor group 
SULi'i) '^SUji{3). As a consequence, the vector and axial currents are classically conserved. How- 
ever, the symmetry observed in nature is S'C/y(3), and only approximatively. This means that the 
classical largest symmetry must somehow be spontaneously broken: SULi'i) SUji{3) -^ SUvi^i)- 
The remaining S'f/y (3) symmetry is also slightly broken due to the different masses of the different 
quarks: 

^^.V^ = + Oimg-nig.), 

df,A^ = + 0{mg) , a / . 

The low-energy spectrum of QCD is made of the well-known octet of Goldstone bosons correspond- 
ing to the eight broken axial symmetries. 

The singlet part of the flavor group deserves a separate discussion, because even in the chiral 
limit, the axial current is not conserved due to the presence of the axial anomaly. Nevertheless, 
the anomalous terms are proportional to the inverse of the number of colors Nc, so Ua{^) can be 
indeed considered as a good symmetry if one also allows Nc to be large: 

a^< = + OK) + o(i-) . 

In this case, one can also think of this symmetry as being spontaneously broken and a ninth 
Goldstone boson appears. 

Following the spirit of Chiral Perturbation Theory, one is led to build an effective theory con- 
taining nine pseudoscalar particles, associated with the spontaneous breaking of the symmetry 
t^L(3) Uii{3) -^ Uy{3) and conveniently collected for this purpose in a 3 x 3 matrix U: 

U{x) = exp i ^ " , where {Xa}a=o,...,8 are the C/(3) generators. 



a=0 



f 



As usual the external sources Sa{x), Pa{x), v^{x) and a^{x), coupled to the vector, axial, scalar 
and pseudoscalar QCD currents respectively, are introduced in order to generate Green's functions 
for the QCD currents and because their behavior with respect to the flavor group transforma- 
tions, taken from QCD, guarantees the reproduction of the QCD symmetry breaking pattern in 
the effective theory. For instance the explicit symmetry breaking effects due to niq ^ are in- 
troduced by freezing the value of the scalar source s{x) = 2BAi, where i? is a constant and 
Ai = diag{mu, ma, rus) is the quark-mass matrix. Similarly, in the U{3) U{3) case, the source 
6{x), coupled to : asG{x)G{x) : in the QCD Lagrangian, provides an excellent tool to keep track 
of the effects of the anomaly. 

In terms of these sources, the divergence of the singlet axial currents reads: 

#^ = 2,/^( E -^«^%^ - «/^) • (3) 

6df,a^{x) V"/^ a ^Pa{x) ^ 69{x) J ^' 

The non-conservation of Aq originates from two different sources: firstly from the quark masses 
through the pseudoscalar currents coupled to Pa{x) and secondly from the anomaly through the 



current coupled to 6{x). When Sqcd is replaced by "S^^pt, (|3|) is automatically satisfied (the sources 
have been defined to do so!). Recall however that, in this case, the equation refers to infinite series 
of the fields, because the different QCD currents are represented by infinite series of pseudoscalar 
mesons. 

These two building blocks (the matrix U and external sources) and the symmetry constraints 
allow for an infinite number of operators in the Lagrangian, but they can be classified depending 
on the associated power of energy (one derivative ~ p ~ ruq ). The complete 0{p^), O(p^) and 
0{p^) Lagrangians were given in R. The problem is that each term in the Lagrangians can be 
multiplied by an (almost) arbitrary function of a special combination of the singlet field and the 
external source 6 given hy X = i ^ 4>o + i 6. Fortunately, the Nc counting allows one to think of 
these functions as infinite series in X whose coefficients would be suppressed in one power of l/Nc 
for each power of X. Due to discrete symmetry restrictions, these series must be either made of 
exclusively even or odd powers of X. Following the notation introduced in [0], Wk{X) (k=0 to 6) 
will refer to the arbitrary functions that multiply 0{p^) and 0{p'^) operators, and Lk{X) (k=0 to 
57) will be used for the ones that multiply 0{p'^) operators. The coefficients for the expansion in 
powers of X are the parameters that are used in practice: 

/ 



W3{X) = -iL(vi,iX + V3,3X^ + 



4 
Lk{X) = Lk^Q + Lk^2X'^ + Lk^^X^ + ..., for even operators, 

Lk{X) = Lk^iX + Lf^^X^ + ..., for odd operators. 

This strengthens the idea that the only possible way of working within this theory is thus 
through the use of two simultaneous expansions in powers of masses and momenta and in powers of 
l/Nc- Luckily simple arguments based on the value of the r]' mass |Q] suggest that both expansions 
can be merged by assuming: p"^ ~ M'^ ~ niq ~ l/Nc ~ ^ (M being the typical meson mass). 

Under this assumption, the leading-order Lagrangian reduces to three terms: 

where: 

D,,U = d^U-ir^U + iUl^, 

X = 2B (s + ip), 

and 

tO 



vo,2 ~ ^ , / ~ ViVc , Br^N)!. 



Notice that the actual order in 5 of any particular computation will depend on the number of 
external fields, because each field carries a 1// ^ Nc ^ (^^' ^ factor. 



A simple dimensional analysis shows that the one-loop diagrams will always be suppressed by a 
factor M^//^ — which is ©(iJ^) according to our choice of 6. As a consequence, the next-to-leading 
order Lagrangian, that is suppressed with respect to the leading one by a factor 6 only, can also 
be treated classically because quantum corrections would only appear in the third order in the 
expansion. 

The next-to-leading Lagrangian involves new 0{p'^) and 0{p^) terms: 

j 

where j labels the 0{p^) operators Oj whose coupling constants Lj are of 0{Nc) (the number of 
required 0{p'^) operators depends on a choice that will be discussed in section 4). Either 114^0) ^5,0 or 
VQfi can be set to zero by means of a change of variables; we shall use v^^ = and v^^i, v^fi, vq^ ~ 5. 



3 Masses and interpolating field in the isospin violating case 

The first step in the calculation is necessarily the identification of the physical states contained 
in the theory, i.e., the diagonalization of the 0{4P') effective action. This was already done in the 
isospin limit in H. However, since the rj/r]' — > tttttt decays are isospin- violating processes, the 
amplitudes for these processes must be evaluated in the case m^ / md / rus- 

In the next-to-leading order, the two-point functions are described by an 0{(jP') and O^S"^) 
effective action given by (^: 

5^2 = - / d^X ( df,4)a Aab d^Cpb - (paBab 4>b) , (5) 

In the isospin-violating case, 

A = I + /\A + I\dA , B = 2Bm{D + dD + AD + AdD) , 

where m = {rriu + md)/2. The suffix A indicates the next-to-leading terms. The terms with the 
suffix d are proportional to nid — rriu- The matrices D, AA and AD Q do not include isospin- 
breaking contributions [|: 



Dn = D22 = D33 -- 


= 1 


L»44 = 


= -D55 = Dqq = D-j-j 


X 


D88 = l + -X , 


^08 = 


V2 
3 " 


, Doo = l + -x- 


3 W02 
2 mB '' 



^Notice that we have chosen a sUghtly different set of parameters than in IS]; now v^fi = and va,o 7^ 0. 



The strange quark mass is introduced through the quantity x = {itIs — m)/m. 

The 0{p^) operators involved in the 2-point functions are O5 and 0%. There are also two 
0{p'^) contributions, because v^^i and ^4^0 are of 0{1/Nc). All these terms contribute to the 0{5'^) 
corrections: 



16 1^5,0 
J 

16 L^fl „ ( . X 



A^n = A^22 = A^33 = —PT^ Bm 

AAu = AA55 = A^66 = A^77 = ^^^jr^ Bmn + ^ 



16 L5,o „ {^^2 
-j^Bm[l + -: 

1 6 L5,o V2 
P 3 



A^88 = — 72^ Bm[l + -x] , 



A^os = Tj^ — Bm X , 



and: 



AAoo = — -^ Bm(l + -x] - 3 ^4,0 , 



ADu = AD22 = AD33 = ^^^ Bm 



AD44 = AD55 = AD66 = AD77 = — ^ Bm(l + x+-x^] , 
ADss = "" .7'" ^"^ ( 1 + ^ a; + ^ x^^ , 



32_L8^o D_ /^i , 4 ,2 2 

AI)o8 = -^^^ ^Smx(2 + x) + V2 z;3,i x , 
f^ 3 

AAo = ^^ i^rri (^1 + ^ x + i x^) - 2 z;3,i (3 + ^) 



When m,u / wid, the following (leading and next-to-leading) pieces must also be considered: 

P 



O J 

Ad^44 = AdA^^ = —AdA^Q = —AdAjj = — -^^ B {rriu — m^) 



Advlos = .2' V 3 "^ ''"^" ~ "^'^'^ ' 



d-D44 = dD^^ = —dDQQ = —dDjj 



rriu -md 



4 m 
1 mu-md 



dD38 = 


V3 2m ' 


dDo3 = 


[2 mu-md 
V3 2m ' 
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AdZ)44 = AdZ)55 = -AdZ)66 = -AdDrr = ^B ^""^ ^'''^ 



2 



J2 ^^2 



AdD3s = — j2^ -y= B [rriu - ma) , 



3 2 Lsfl /2 
A simple change of variables: 



AdDos = — j^ \l7;B {niu - nid) 



V = (/ + ^A^ + iAdA),/. (6) 

provides the correct prefactor in the kinetic term and reduces the calculation of the masses to an 
eigenvalue problem. The matrix to be diagonalized is: 

2Bm(D + AD- -{AA, D} + dD + AdD - -{AdA, D} - -{AA, dD}) . (7) 

This matrix is diagonal in the kaon and the charged-pion sector, so the interpolating fields and the 
particle masses can be written straightforwardly as: 

1 1 

7r+ = ^(V'l + #2) , 7T- = -^{tpi - #2) , 

K+ = — y={tp4 + ii^b) , K- = —=0)4 - iip^) , 
Kq = — ^(V'6 + #7), Ko = y={^e - it/jj) . 

Ml^ =Ml_=B (m„ + mrf) (1 + 8 5 (m„ + m^) ^ ^''° ~ ^''° 
Ml+ = Ml- =B{mu + nis) (1 + 8 B (m„ + m,) ^ ^^'° ~ ^^'° 
^io = ^h = B{md + ms) (1 + 8 S (md + m,) ^^^'° ~ ^''° 

ttq, r/ and rj' are related to (p^^ (po and (j)8 by a rotation S*. In the rriu = rrid case, the rotation 
mixes r/ and ry' only: 

/I 

5=0 COS0 sin< 
V — sin cos ( 

The isospin-violating terms will be considered as small corrections to be treated perturbatively. To 
first order in this expansion, the relation between the original fields and the eigenvectors is given 
by a matrix F ^■. 

VJ> = iF^)pa(pa, With F<' = {S'')-\l + ^AA + ^AdA), 



t 



These expressions are to be understood in the sense of perturbation theory: 

F'' = S''' (I + ^AA + ^AdA) + (IT S''' (I + ^^A) + AdT S'' + 0(5'^ ,(f 



and S"^ = S{I-dT- AdT) . 
The matrices dT and dDT are defined by: 

dTpQ = ^2 _j^2 i "PP I dD\ ifg), (v? = S^;) , 

^dTpQ = -lE^{^p\AdD-^{AdA,D}-^{AA,dD}\^Q). 

The only non-vanishing elements in these matrices are: 

dTT^n = —dTjjTT = eo(cos^ — v2sin0) , 

A r- 

dT^ri' = -dTr^i-n = 2 _ f^2 ^0 (^™ ^ + V 2 COS 6) , 

AdT^^ = -AdTriT, = (e3,8cos6'-eo,3sin6'), 
AdT^n' = -AdTn'T, = ^ _ j^f 2 (ea.ssing + eo^scosg), (8) 



where 



€3,8 = eo 16 (M^ - Mk) -2 , 

eo,3 = eo V2 (^ - 3^3,1 + 12^0,2 ^ + 16 (M; -M^) '—^ ^ + 2 ^4,o 



2 I --- V-'^TT '-/i. / (•2 



A = M^-M^ . 



The dimensionless quantity eo is a good measure of the isospin-breaking perturbation since 
B{md - niu) = -eo A VS (1 - 16 M^ 



2-Z>8,0 — ^5.0 



and can be estimated in terms of the observable quantities A and M^ (^): 



eo = -^, Ml = {Ml^-Ml^)-{Ml^-Ml^). (9) 

The combination of masses Ml isolates the QCD isospin-breaking effect due to the quark masses, 
because the electromagnetic contribution to the mass splitting is the same for the kaons and the 
pions and will cancel out. 

It can be checked that the leading-order contributions to the matrix F do indeed match the 
diagonalization given in p. 

As a consequence of (^ being the only non- vanishing terms the first isospin- violating corrections 
to M^o; ^v ^'^d -^n' ^^^6 of ^(^0) which goes beyond our working precision and will be therefore 
neglected. 

In any case, for the isospin-violating decays that are the subject of this paper, the amplitude is 
proportional to eo so only the complete isospin-violating eigenstates and the complete quark mass 
matrix are required for the calculation. The free parameters in the theory can be estimated in 
the isospin limit since any eo correction would give a second-order contribution of ©(eo) to the 
amplitude. 



4 Four-point processes in Ul{3) ^i?(3) xPT 



To leading order, the only two relevant four-field terms in the Lagrangian are formally the same as 
those that appear in the SU{3) theory: 

■ fabr fcdr 4>a d^(t)b (t>c d^4>d + l^TJo'^o-br ^res dgcd 4>a (t>b (t>c4>d'^B Me ■ 



The next-to- leading 0{S'^) contributions can originate from: 

a) terms with constants of 0{1/Nc) from the 0{p^) Lagrangian; 

b) terms with constants of 0{N^) from the 0{p'^) Lagrangian; 

c) terms with constants of 0{Ne) from the 0{jp'^) Lagrangian. 

The A^c-power counting for the coupling constants in the model was studied in Q- According 
to this work there can be no 0{jfi) correction, and there will be only one more 0{p'^) contribution, 
associated with the coupling constant v^i. The next chiral order looks less encouraging: at first 
sight one might think that nine independent 0{p^) terms have to be included (Oi, O2, O3, O5, Os, 
Ois, Oi4, O15 and Oig). This problem can be dodged by noticing that the constants L, associated 
to most of them are 0{Nc) due to the contribution of the term Oq that was eliminated through 
the Cayley-Hamilton theorem. A more convenient set of independent operators can be chosen by 
eliminating Oiq instead. In this case, the only operators left in the list are Oq, O3, O5 and Og, 
because all the other terms are suppressed by a factor 1/Nc or more. 

The new set of coupling constants {Mi} is related to the old one {Li} . In particular, for the 
first three of them, one obtains : 

Mifl = Li^o-^ ; M2,o = L2,o - Mofi ; Mg^o = ^3,0 + 2Mo,o • 

One can use the first two relations to get an estimate of the new constant Mq^q. Mi^q and M2,o are 
expected to be both of 0{N^), and thus negligible in front of Mq.q, -^1,0 and ^2,0- Furthermore, 
as shown in [10|, Li^, L2,o and L3.0 are equal to the corresponding Li, L2 and L3 from the octet 



theory up to one-loop corrections — which is beyond our working precision. Ls^q will be therefore 
directly borrowed from the 5C/(3) model (see, for instance, |11| and references therein): 



L3,o = (-3.5 ±1.1) -10-3. (10) 

Li and L2 can be used to fix the new constant: 

Mo,o = ^(^^1,0 + ^2,0) + 0(iVc ) - (1-2 ± 0.4) 10-3 . 



Another way to estimate Mq^o is given by the QCD bosonization models |12]: 



'>*'» = il^-^i-^**-'""' 
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This value seems more reliable since bosonization models have proved to give excellent results for 
the 0{p^) constants (except for the operators that contain explicit symmetry breaking because 
then the results are model-dependent, but this is not the case for Oq). 

Yet a third possible source to fix Mq^o can be used: it turns out that the most important 
contributions to the decay ry' — > r/vrvr come from Oq, so the experimental data for this process can 
be used to determine Mq^ (see section 7 for details). The quoted error is the one induced by the 
experimental uncertainties: 

Mo,o = (1.54 ±0.1) -10"^ . (11) 

This value is in good agreement with the theoretical estimations discussed above. It is worth 
pointing out that M^^ turns out to be quite small. Using the central value in (0), we take: 

M3,o = -0.4 • 10-3 . (12) 

The other constants that appear in the calculation are B m, x f\ , f, L^^q, L^fi, v^^i, vo,2 and 
t;4^o- Their values are fixed by the masses and decay constants Mj^, Mk, M^, M^/, /,ri Ik, frj and 
frj'. Within the required precision and in terms of the mixing angle 9, the first correction can be 
written in terms of measurable quantities by using the following identities: 

8 , n n, , , M2 + 3M2 - 4M|. + 3(M2 - M2) sin^ 9 
AM = ^(Mi-M^)(2L,o-L,o) = ^ ^ 4(Mi-M^) ' 

12 3 {M^,-M^) 

^N = 3 .3,1 - J, vo,, L5,o = 1 + ^ jMl^Mf) ' 

-^5,0 ^ 1 / Ik _ ,\ 

/2 4(Mi-M2)lA J' 

..o^-^(^-2). (13) 

The parameters from the leading-order Lagrangian must be evaluated up to 0{5) corrections: 

X = 2j^{1-Am)-2{1 + Am), 
f = /.(l-^M^^), 

-3t;o2 = (m2 - ^^^^i±i^) (1-3^40) + ^ (Mi -M2)(1 + Am -A;v-^.4,o) tan 
- ^((Mi - m2) Am + {2Ml + M^) An) . (14) 

As shown in [^, the fitting does not fix the value of the mixing angle 9. One observes, however, 
that the corrections on M^ and M^, given by Am and Atv are minimized for 9 between —20° and 

^Recall that, to first order in m„ — m^, the isospin-breaking eflect is an overaU factor in the ampfitude. 
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—22°. Furthermore, this minimum sensitivity prediction agrees with the experimental data and is 
reasonably close to the leading-order prediction (9 ss —21.7°), as expected if the U{3) expansion is 
to make sense. 



5 The T] -^ TTTTTT decays 

A first estimation of the decay rates can be given by the current algebra |1S] or by the leading order 



0{p ) in SU{3) Chiral Perturbation Theory [^], Q. According to this theory, the amplitude is: 

Ao+-{s,t,u) = -^is-^M^). 

Electromagnetic contributions to the decay amplitude need not be included to leading and next- 
to-leading order in the low-energy expansion, as argued in [^]. The only relevant QED effect in 
this case is the splitting between the masses of the charged particles, which is taken into account 
through the definition of eo (w • 



The numerical results 


are 


too small: 














-^ VroTToVTo) = 

> 7ro7r+7r_) = 
r = 


= 100 eV 
= 66 eV , 
= 1.51 . 



The predicted rates are off by a factor of four. The branching ratio r is however not that bad. 
This is not such an amazing feature, since it has been shown that the ratio is related to the isospin 
selection rules (see comment on section 1 and appendix A). In any case, it seems to indicate that 
whatever is missing in both calculations is essentially a multiplicative correction that cancels out 
in the ratio. 

The next-to-leading 0{p^) calculation in SU{3) Chiral Perturbation Theory was carefully ana- 
lyzed and presented in 0: 

r'5f^(3)(^^7ro7ro7ro) = 229 eV , 
r^f^(3)(^^7ro^+7r_) = 160 eV, 
r = 1.43 . 

The new terms come from the 0{p^) Lagrangian and from one-loop diagrams. An unexpectedly 
large contribution arises from the unitarity correction term, especially from the pieces that corre- 
spond to vrvr final-state interaction []. Technically, these contributions turn out to be so important 
because the one-loop diagram integrals are plagued with infrared divergences that produce a sig- 
nificant enhancement of the quark-mass perturbation. 



^The importance of the so-called rescattering effects had been already pointed out in some previous work made in 
the context of current algebra which attempted to improve the initial result by including the tttt final-state interactions 
through the imposition of unitarity and analyticity |15[| . 
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5.1 Leading order in UiiS) Ur{3) xPT 



As shown in section 3, the leading-order Lagrangian required to describe this process is very similar 
to the Lagrangian used in SU{3) xPT. The difference with respect to that case lies in the mixing 
effects. The anomaly appears only indirectly through the contribution of fo2 to M^ and 9. 



An 






({cose-V2sme){s--M^' 



M2 r- 

H tan 6* (sin 6* + 4 V2 COS 6* 



(15) 



The last term is proportional to M^ so it produces a small correction to the final result. The mixing 
effects reduce essentially to a factor of cos9 — \/2sm9 (=1.45 for 9 = —21.7°) ^, 16|. 
The predicted rates are higher than in the octet model: 



r(r?- 


-^ vrovrovro) 


= 180 eV, 


r(r?- 


> 7ro7r+7r_) 


= 121 eV, 




r 


= 1.49 . 



5.2 Next-to-leading order in [7^(3) ® Ur{3) xPT 

The final result for the squared amplitude (consistently expanded into leading-order piece -|- cor- 
rections) is a very lengthy expression that can be found in the appendix. In particular, most of 
the free parameters in the theory can be expressed in terms of physical quantities (13, |lj). Only 
Ls^o and Mofi have to be estimated numerically ( p!0| , 11). The final expression is a function of the 
mixing angle 9, because the determination of the free parameters through the masses and decay 
constants does not fix the value 9, although there are many reasons (see end of section 4) to expect 
it to be around —21°: 



9 


-20° 


-21° 


-22° 


r{r, 


-^ TToVroTTo) 


129.1 eV 


124.8 eV 


119.7 eV 


nv- 


-^ VroTT+TT-) 


87.9 eV 


84.7 eV 


80.9 eV 


r 


1.47 


1.47 


1.48 



The predicted rates have improved but remain too small when compared to the experimental 
data. This result seems to point out that the pion-pion final-state scattering effects are indeed 
a key element in the process. In the context of C/l(3) (8) Ur{3) xPT, these are next-to- next-to- 
leading order contributions. When properly evaluated, this should include the one- loop correction 
that stems from the leading-order Lagrangian, and a myriad of new terms — including some 0{p^) 
terms — coupled to totally unknown parameters. One expects, however, that the most important 
contribution should come from the pion-pion interaction and that it should not be much different 
from the SU{3) result @. 
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6 The rj -^ tttttt decays. 

The calculation goes along the same lines as the one presented in the previous section although 
there is obviously no SU{3) prediction to compare with. The processes are also isospin violating 
and will be evaluated to first order in eo- 

The expressions (ll]), (0) and the vector-meson dominance argument also apply to this process, 
so the ratio between the charged and neutral channels can also be predicted to be around 1.5. This 
number is compatible with the experimental data. 



6.1 Leading order 



The leading-order amplitude for the rj' -^ 7ro7r_|-7r_ decay reads: 

4 



An 



en „ sin 6* + \/2 COS /, ^ , ,9, , ^ 

" tan 6* ^ ({s--M^) {V2 cos 9 + sinl 

V2cos0 ^ 3 



sin( 



— ^ cote (4\/2 sine- COS e 



(16) 



Notice that 6^9 + ^ exchanges rj and rj' (except for a global sign) only in the isospin limit 
because the it — tj and vr — r/' mixings are not equal, so this symmetry cannot be used to relate (^ 
and (IC). 

The integration over the allowed phase-space region leads to: 



nv' 


-^ TTq-Kq-Kq) 


= 457 eV 


V- 


■> 7ro7r+7r_) 


= 405 eV 




r 


= 1.13. 



The prediction is of the correct order of magnitude. These results agree with the estimations in 
[17|. Notice, in particular, that the present computation does include the effects of tjq — r/g mixing 
and the gluon anomaly that they point out as a crucial ingredient of their calculation. 



6.2 Next-to-leading order 

The analytical expressions for the amplitudes are given in the appendix. It turns out that the 
numerical results are not good, because the corrections are huge: 



9 


-20° 


-22° 


T{r]' -^ vroTTovro) 


2280 eV 


2137 eV 


T{r]' -^ 7ro7r+7r„) 


1642 eV 


1536 eV 


r 


1.39 


1.39 



One would a priori expect close similarities between both rj and rj' — > tttttt processes. The rj' is 
however nearly twice as massive as the 77, so the emerging pions will have relatively high momenta. 
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This might suggest that the bad results quoted above are to be blamed on a breakdown of the 
low-energy expansion. However, high momenta also allow for a lot of phase space for pion-pion 
interactions that produce intermediate resonances like the p or the a. This suggests that, in this 
case, the effect of resonances might be even more important that in the r/ — > vrvrvr. This would 
correspond to one-loop corrections, but one can barely expect a third contribution to cancel the 
huge second-order corrections in order to produce a good final result. 



7 The Tj' -^ rjTTTT decays. 

This transition is also described by the Lagrangian discussed in section 4 and will not include new 
elements to be taken into account. As a matter of fact, the calculation is simpler because these are 
not isospin-violating processes so one can set eo = everywhere. 



7.1 Leading order 

The leading-order amplitudes for both the charged and the neutral channel read: 

A^.^^,Js,t,u) = ^(2V2cos(20)-sin(20)) . (17) 



The decay rates that follow from ( [T7| ) are however very small [16, |^, 19|: 



r(r?'- 


-^ r/vrovro) = 


= 1.0 keV 


r(r?'- 


>777r+7r_) = 


= 1.9 keV 




r = 


= 1.9. 



The actual values for the rates are however 40 times smaller than the experimental ones. A 
possible justification for the leading-order prediction being so small was already pointed out in 
[19, ^: notice that the leading-order amplitude (17) vanishes in the chiral limit. 



7.2 Next-to-leading order 

The first non-vanishing contribution in the chiral limit would come from the next-to-leading La- 
grangian that does not break chiral symmetry explicitly: Oq and O3. The actual computation 
shows that even when massive quarks are considered, the most important contributions to the am- 
plitude do indeed come from these operators. Numerically the contribution from O3 is suppressed 



by the small value of M^^ (12). In the center of the Dalitz plot, where s = t = u 



M^,+M3,+2 Ml 



17' 



V 



3 



the contribution from Oq to the total amplitude is equal to 10 times the leading-order amplitude. 
This seem to indicate that the key mechanism associated to this transition is related to some qual- 
itatively special dynamics that are not reflected in the leading-order Lagrangian, but show up in 
the next order (and, in particular, in the Oq operator). As in the previous sections, the analytical 
expression for the amplitudes has been relegated to the appendix. 
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As a consequence of this strong dependence on one particular operator, these processes provide 
an excellent way to fix the unknown constant Mo,o- In this work, the fitting was done in the charged 



channel (11) 



From a strict perturbative point of view, only linear next-to-leading corrections should appear 
in the decay rate. This implies that the squared absolute value of the amplitude is introduced in 
an expanded form: |AioP + 2|^lo||^A'"Lo|- In the case that is considered in this section, however, 
this expansion does not make sense, since the second-order contribution to the amplitude is by no 
means a small correction to the leading-order value! The relevant dynamics appear in the next-to- 
leading Lagrangian and should not be treated as a perturbation. The decay rate must be computed 
from the square of the absolute value of the amplitude, \Aio + ^nloI"^, and not from an expanded 
form. Otherwise one obtains a negative value for the rate! This is not a sign of a poorly convergent 
expansion, because third- and higher-order corrections could still be expected to be small. The 
point is that the leading-order contribution is (nearly) zero, so the next-to-leading contribution is 
actually to be considered as the first term in the expansion. 

The results are also a function of the mixing angle, but the dependence is negligible in the 
region of interest, -20° <9< -22: 



r(^'- 


-^ r/TToVTo) = 


= 48.7 eV 


V- 


^ ?77r+7r_) = 


= 88.9 eV 




r = 


= 1.83 . 



The agreement between the predicted and the experimental value for the neutral channel is 
almost within the experimental error. The ratio between the two channels also matches the mea- 
sured value. A word of caution is however needed in relation to these results: the decay rates turn 
out to be extremely sensitive to the value of this constant. Some illustrative figures: a deviation of 
+0.05 • 10~^ in Mqa increases the estimated rates by 30%. The ratio remains almost unaffected, as 
dictated by isospin symmetry. 



Some recent work on a pseudoscalar-scalar meson coupling model [^] (see also [22|) indicates 
that the dominant contribution to rj' -^ rjinr decays comes from the exchange of the intermediate 
scalar resonance ao(980). It also predicts a less important contribution from the a, and no tree-level 
contribution from the p (forbidden by G-parity). These features agree from a qualitative point of 
view with the results in this paper, since the ag resonance does indeed contribute to Mq^ [p^ j. 
The p internal exchanges are a delicate issue in the f7(3) formalism, because M^/ > Mp, which 
means that the p particle has not been fully integrated out. The G-parity constraint protects the 
r]' — > r]Tnr system against this obstacle. The only relevant (but smaller) one-loop corrections should 
correspond to a exchanges. 

8 Conclusions 

The Ul{3) Ur{3) formalism provides a systematic way of dealing with the r] — ri' — ir interactions. 
In particular, it is supposed to give a fairly accurate description of the rj/rj' system. The good 
results obtained in [^ for the values of the masses and the mixing angle 9 are also certainly a 
strong motivation for the study of these decays. 

The effective theory is built upon the assumption that all particles other than the nine Gold- 
stone bosons can be integrated out, so their effects would show up in the effective vertices. This 
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assumption is exact in the niq — > and Nc ^ oo limit, because then the bosons are really massless 
and any massive resonance will decouple. As we move away from this ideal situation, quark-mass 
and 1/Nc perturbations are introduced. This produces three sort of corrections. Firstly, new inter- 
actions between the Goldstone bosons themselves appear (consider, for instance, the leading-order 
mass terms, that are either of 0{mq) or 0{\/Nc)). Secondly, the integration of all intermediate 
resonances becomes less reliable. Light resonances like p, to and a can indeed produce serious 
problems, specially when they happen to play an important role in some particular process. When- 
ever this occurs, one-loop contributions, where these resonances appear in terms of interacting 
Goldstone bosons, ought to be essential. This seems to be certainly the case for rj -^ Trvrvr decays 
and probably for tj' — > Trvrvr decays, too. (Notice that the integration of p and lo is more delicate in 
the U{3) formalism than it was in the octet theory, because M^/ is around 1 GeV, which is higher 
than the masses of these particles, so their low modes have not been integrated out. This might 
imply a change in the estimated value of some constants like L^^ because L3 feeds precisely on 
isovector- meson contributions |^3|). In the third place, the anomalous diagrams involving g"^ GG 
— the gluonic part of tjq (see figure) must be associated to next-to-leading contributions, because 
they are suppressed in l/N^ with respect to the diagrams involving quark currents. 



qYsq 





o(N,) ai) 

Figure 1: Quark-current and glueball contributions in r/ — > vrvrvr. A similar diagram can be drawn 
for rj' decays. 

The U{2>) formalism offers a more complete description of the r] than the octet theory, because 
the latter identifies r/g and rj, neglecting mixing effects that should to be rather important because 
the mixing angle is not that small. The f/(3) calculation improves the tree-level SU{3) prediction 
for the ry — > vrvrvr decays. The most relevant new feature for the these decays seems to be the 
rj — r]' mixing (the anomaly does play a role there but it is not a direct contribution to the 4-point 
function). The mixing angle is well-predicted at leading order and is stable under next-to- leading 
corrections {9 ~ —20°/ — 22°), so significant corrections due to mixing effects are not expected to 
appear at higher orders. Nevertheless, since the next-to-leading predicted rates stay well below 
the experimental values, there seems to be no way to avoid the importance of unitary corrections 
and low-energy resonances whose effects are not included in the effective vertices precisely due to 
their small mass. The C'((5'^) effective action would include these final-state-interaction corrections, 
but tadpole corrections as well as many new terms in the Lagrangian with unknown coupling 
constants should also be taken into account. If the discrepancies with experimental values are 
indeed mainly due to the presence of intermediate states, tadpoles and 0{6^) counterterms could be 
safely neglected and the unitary corrections — that can be computed with the known ingredients — 
should be the only contribution that matters. This one-loop calculation is however out of the scope 
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of the present article and will be left for future work. (A first step in this direction is the evaluation 



of the one-loop contributions to the masses and decay constants |24]). 

The same discussion should also apply to the r/' — > vrvrvr decays. Although the leading-order 
approximation produces good estimates of the rates, the emerging pions are far from being soft and 
this is probably the reason why the expansion apparently blows up when corrections are included. 

rl — > r/TTvr transitions are probably the most interesting decays in this paper, because this decay 
could not be analyzed in the framework of S'C/(3). Momenta are not expected to be too high, so 
the low-energy expansion might actually work. The results are reasonably good at tree level, but 
an estimate of the one-loop corrections would also be of great interest in this case, in order to check 
the convergence of the expansion, which is certainly one the most fragile points in the ?7(3) theory. 
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A Isospin symmetry for 77/77' -^ tttttt decays 

The relation ([l|) between the charged and the neutral channels can be proved by simply considering 
isospin symmetry. By definition, for the charged channel, 

A{s,t,u) = (vro(pi) vr+(p2) 7r_(p3) I T/(p4)) ; 
A{t,u,s) = (vro(p2) vr+(p3) 7r„(pi) I ??(p4)) ; 
A{u,s,t) = (vro(p3) vr+(pi) 7r_(p2) I ??(P4)) • 

The permutation of momenta can be viewed as a permutation in the isospin values instead. For 
instance ||: 

kofe) TT+iPs) vr_(pi) ) = -\pi p2 P3 ; 11 1-1 10 ) ; 
Kofe) vr+(pi) 7r_(p2) ) = -\pi P2 P3 ; 11 10 1-1 ) ; 



These three particle states can be expressed in the Clebsch-Gordan basis |/(^'^) / M) where the 
states are labeled in terms of the isospin from particles 2 and 3 I^'^''^> , the total isospin / and the 
component in the z-direction of the total isospin J^. 

"The minus signs are due the fact that (71"+)* — — tt-. 
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Due to the bosonic nature of pions, only the totally symmetric states will contribute to the 
decay amplitude. This restricts the final state to a superposition of three states (expressed in the 
Clebsch-Gordan basis): 

6 4 2 

\symmetric state) = ^=|2 3 0) + ^=|2 1 0) + -^\0 1 0) . (18) 

The neutral-channel analysis is much simpler, since all particles are / = 1, /^ = 0, which can 
only produce totally symmetric states: 

110 10 10) = y||2 3 0)--|=|2 10)--||0 10). (19) 

However, the rj -^ tttttt transition must be induced by the isospin-violating piece in the QCD 
Lagrangian |||]: 

>Cqcd = -^C"^" - "T-rf) ("""" - dd) . 

This is a A/=l operator, so the 1=3 pieces in ([iq) and (|l^) will not contribute to the decay 
amplitude. The remaining terms differ by the value of /(^''^) , so this number can be used to label 
the two different contributions to the charged amplitude: 

A{s,t,u) = ( 10 11 1-1 I 00 ) = -^ ^0 - -^ ^2 ; 

A{t,u,s) = { 11 1-1 10 I 00 ) = y|- A2 ; 

A{s,u,t) = A{s,t,u) ; 

A{u,t,s) = A{u,s,t) = A{t,s,u) = A{t,u,s); (20) 

and to the neutral amplitude: 

A{s,t,u) = (10 10 10 I 00) = --^Ao--^ ^2. (21) 



From ( poD and (pTD, it is straightforward to check that: 

A{s,t,u) = A{s,t,u) + A{t,u,s) + A{u,s,t) . 
Obviously the analysis applies to both rj and r/' decays. 

B Next-to-leading decay amplitudes 

The analytical expressions of the decay amplitudes are lengthy due to the dependences on the 
mixing angle 0. As discussed in the main body of this article, every constant in the computation 
except Mq^o ^-nd L^^q can be written in terms of measurable quantities and 9: 

3 M2 + M2 - 4 Ml- + 3 sin2 (M^, - M^) 
^'' = 4 (Mi - M2) ' 
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3 M^, - M2 
Ay/2 Ml- M^ 

Ap = (^-1); 

V,, = -2^ + A' o. 



(22) 



The amplitudes will always have the general form: 



ANLo{e) = ALo{e) +^M{0)AM{e) + ^N{e)AM{9) + ^pAp{9) + vi^QAi,^{e) 

The actual expressions for these next-to-leading contributions to the charged-channel decay 
amplitudes are given below. The amplitudes for the neutral channels can be inferred from them. 

B.l 77 — !► TTTTTT 



A 



M 



An 



3cos2 6' (^/2 cos 6' + sin 61) 2 (cos 61- ^/2sin6') 

3s cos^ 6* (-2 + 7 cos^ 6-7 cos"^ 6* - 2\/2 sin 6^ cos 6' + 4 sin 6^ cos^ 6) 
+ Ml cos 9 (35 cos 9 - 104 cos^ 6^ + 85 cos^ 9 + 5\/2 sin 6^ + 12\/2 sin 6^ cos^ i 

- 37^2 sin 6' cos^ 61) + -—^^^ (-2 - 40cos2 6* + 106cos'^6' - 64cos^ 6* 

M^-M^ 

- 13\/2 sin 6' cos 61 + 4\/2 sin 9 cos^ 6* + 25\/2 sin 6^ cos^ 9 

1 



3(^/2 cos 6* + sin 6*) (cos 6* - ^/2 sin 
3s sin 26* (cos 26^ + 
+ 10M2sin6' (3-5cos2 



Aj 



3s sin 26* (cos26l + ^ sin26') + M^\/2 (l + 16cos^6l- 17cos'^6') 



Ml 



9cos2 6'sin6l (sin6' + \/2cos6')2(cos6l - ^/2sin6l) M^ - M^ 
3s (-2\/2cos6l + 5\/2cos^6l-2\/2cos^6'-\/2cos'^6l 

- 8 cos^ 6* sin 6^ + 26 cos^ ^ sin 6* - 22 cos^ 9 sin 9) 

+ Ml (21\/2 cos 9 - 37\/2 cos^ 9 - 13\/2 cos^ 6* + 29\/2 cos^ 6* 

+ 2 sin 6* + 72 cos^ 6* sin 6* - 210 cos"^ 6* sin 6* + 152 cos^ 9 sin 9) 

Ao,o = 8 sin 9 (4\/2cos 9 + sin 9) fs^ - 3sso - ut + 3M| 
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^3,0 = -, ^=— — ( cos6' (s -15sso + 2tu) 

3 cos t^(cos y — V 2 sm 6) ^ 

- cos6lsm2 6l (42M| + 13s^ + 57sso - 16tn) - 2\/2sm6l (3M| + s^ + 3sso - tu) 
+ 2\/2 sin e cos^ 6* (15M| + 4s^ + 30sso - 7iu)) , 

where so = M^/'i + M^ and M^ = M^ (M^ + M^). 

B.2 ?]' ^ TTTTTT 



A 



M 



9 COS 61 sin 61 (cos 61- ^/2 sin 0)2 
9s (-3 cos 6* + 8 cos^ 61-5 cos^ 6 - \/2 sin6l + 2\/2 cos^ 6* sin - \/2 cos"^ 6* sin ( 



4 M'^ 
+ —2 ^ (-16 cos 6* + 56 cos^ 0-38 cos^ 6 + 2\/2 cos^ 6* sin 6* - 13\/2 cos"^ 6* sin 6*) 

+ M^ (92 cos 6* - 267cos^6' + ITScos^^ + 28^/2 sin - 71\/2cos^ 6'sin6' 
+ 62v/2cos''6'sin( 



An = p (mI (\/2 cos^ 6* - 10\/2 sin^ 61-28 sin 9 cos 6*) 

3(cos6' — v2sin0) ^ 



Ap 



+ 3s sin 9 ( V2 sin 6* + 2 cos / 

4 M^ 



9 sin 9 cos 6l(cos - ^/2 sin 9f M^, - M^ 
'mI (- 16 + 108 cos^ 6 - 246 cos"^ 9 + 152 cos^ 9 - 24\/2 cos 9 sin 

+ 74\/2 cos^ 6* sin - 29\/2 cos^ sin 9) - 3s(-4 + 22 cos^ 61-40 cos"^ + 22 cos^ 9 
- 2\/2 cos sin + 5\/2 cos^ sin - \/2 cos^ sin 0)) , 

^4,0 = --^ An , 

8(sin0 + \/2cos0)(4\/2sin0-cos0) / 2 2\ 

Ao,o = ^ ^ '^^ . ^ '- (s^ - 3sso -ut + 3M| , 

COS0 — v2sin0 ^ ^ 

^3,0 = -, 7=— -7 ( - Sin (s2 - 15SS0 + 2tu) 

3cos0(cos0 — v2sin0) ^ 

+ sin0cos^ (42M| + I3s^ + 57sso - 16tn) - 2\/2sin0 (3M| + s^ + 3sso - tn) 

+ 2\/2 cos sin^ (15M| + is'^ + 30sso - 7tu)] 

where sq = M^,/3 + M^ and M| = M^ (M^, + M^). 
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B.3 r^ -^ rjTnv 



i-M 



9cos0 M^-M^ 



(2\/2 cos^ 6* + 2 sin^ 6* - 3\/2 cos 61 sin^ (9) , 



M2 r- 

An = — ^ (V2cos26l-2sin26l) , 

2 M^ / 

9cose(^/2cos0 + sin0) M^ - M^ \ ^ ^ '? '' ^^ 

+ 4 cos^ e (3M2 - SM^, - 4M2) + ^/2 cos sin (3Af 2, + 7 Ml) 
+ \/2 cos^ e sin e (-3M2 + 3M2, - MM^)) , 

^0,0 = -4 (2\/2 cos 26* - sin 26) {m^ + 2M^(M^ + M^,) + M^M^, + s^ - tn + Sssq) , 
. 1 . 

^3,0 — o ^0,0 ) 

where sq = (M^, + M^ + 2M2)/3. 
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